Attosecond neutron Compton scattering from protons 
entangled with adjacent electrons 



C. Aris Chatzidimitriou-Dreismann^ 



^ Institute of Chemistry, Sekr C2, Technical University Berlin, D- 10623 Berlin, 
Germany. Email: dreismann@chem.tu-berlin.de 



Abstract 

The effect of "anomalous" scattering of neutrons and electrons from protons 
in the electron-volt energy-transfer range is considered, and related experimental 
results are mentioned. A recent independent confirmation of this effect with a new 
data analysis procedure is presented. Due to the very short characteristic scattering 
time, there is no well defined separation of time scales of electronic and protonic 
motions. An outline of a proposed theoretical interpretation is presented, which 
is based on the fact that scattering protons represent open quantum systems, thus 
being subject to decoherence. 
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1 Introduction 

Several neutron Compton scattering (NCS) experiments on liquid and solid sam- 
ples containing protons or deuterons show a striking anomaly, which is a shortfall 
in the intensity of epithermal neutrons scattered by the protons and deuterons. E.g., 
neutrons colliding with water for just 100-500 attoseconds (1 as = 10~^* s) will see 
a ratio of hydrogen to oxygen of roughly 1.5 to 1, instead of 2 to 1 corresponding 
to the chemical formula H2O; cf. (T]|2|. The experiments were done at the ISIS 
neutron spallation facility, Rutherford Appleton Laboratory, UK. Due to the large 
energy and momentum transfers applied, the duration of a neutron-proton scat- 
tering event is a fraction of a femtosecond which is extremely short compared to 
condensed-matter relaxation times. 

This new effect has been confirmed using an independent method, electron- 
proton Compton scattering (ECS), at the Australian National University. (3) 14). 
ECS experiments from a solid polymer showed the exact same shortfall in scat- 
tered electrons (with energy about 20-35 keV; scattering angle 45°) from hydrogen 
nuclei, comparable to the shortfall of scattered neutrons in accompanying NCS ex- 
periments on the same polymer. The similarity of the results is striking because 
the two projectiles interact with protons via fundamentally different forces - elec- 
tromagnetic and strong tSlEllll- 
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Due to its novelty and far-reaching consequences, however, this effect has been 
the focus of various criticisms, cf. | 5, 6]. Therefore, considerable work to identify 
possible sources of experimental and data-analysis errors was made during the last 
five years, which succeeded to demonstrate the excellent working conditions of 
the spectrometer Vesuvio at ISIS |7|. Moreover, extending these investigations, 
the complete "exact formalism" of data analysis (3) was applied to NCS-data by 
Senesi et al. [8 |, for the first time; analysis of time-of-flight (TOF) spectra from 
solid HCl revealed the existence of a strong "anomalous" decrease of the scattering 
intensity from H (up to 34%). Additionally, this result was found to be in excellent 
agreement with the corresponding outcome of the standard data analysis procedure 
applied at ISIS |8j. 

Recently, scattering of neutrons in the 24-150 keV incident energy range from 
H2O relative to that of D2O was investigated |9|. In clear contrast to the NCS and 
ECS results, it was claimed that the measured neutron scattering intensity ratios 
exhibit no anomalous behaviour. However, an improved analysis [lOJ of the keV- 
data within the frame of standard theory showed that the considered scattering 
anomaly is present at both 5-100 eV and the keV ranges of neutron energies. 

Very recently, the mentioned standard NCS-data analysis method (7) was suc- 
cessfully compared with a newly proposed (by B. Dorner, ILL) model-free data- 
analysis procedure, the latter being independent of the form of the momentum 
distribution and the resolution function llll . In this work, the original results from 
the metallic hydride NbHo.s f 121 were analyzed. The comparison of results ob- 
tained with the mentioned two independent methods underline the importance of 
the effect under consideration; see Fig. I. 
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Figure 1: "Anomalous" scattering from the solid metallic hydride NbHo.s il2J . Shown 
is the measured ratio Rexp — o'H/o'Nb of scattering cross-sections of H and Nb nor- 
malized with their expected (tabulated) ratio Rconv Broken line: conventional ex- 
pectation; "full" symbols: results taken from the original publications [12J; "open" 
symbols: results of the model-independent data analysis procedure invented by Dorner 
lini(a). 
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2 On scattering time 



In the context of NCS, as provided by the Vesuvio setup, the Impulse Approx- 
imation (lA) is valid 1141 1151 and the characteristic time scale — often termed 
"scattering time", Tsc — of the neutron-proton scattering process is very short 
(3l|71[l2l[ll, 

Tsc ~ 100 - 1000 as (1) 

(as: attosecond). This is a consequence of the large energy (up to 100 eV) and 
momentum transfers attained with the Vesuvio instrument [7| and follows from 
the theoretical result valid in the I A I14II15I 

Tscqvo ~ 1 , (2) 

where vo is the root-mean-square (rms) velocity of the nucleus and hq is the (ab- 
solute value of) momentum transfer from the neutron to the proton. The time r^c 
is given by the i-width of the intermediate correlation function F{q, t), which is 
related to the dynamic structure factor S{q, u) by Fourier transform 

1 

Siq, <^) = ^J exp {-iLjt)F{q, t) dt . (3) 

It is interesting to note that the "actual duration" Tact of a neutron-proton inter- 
action should be even shorter, as a classical estimate indicates. E.g. a neutron with 
kinetic energy Eq ~ 10 eV will pass a distance of 10^^ A (i.e. the range of the 
strong interaction) in a much shorter time. Tact ~ 10^^^ s. In the light of Relativ- 
ity Theory, this has a crucial consequence: the "actual" (or "effective") scattering 
system — i.e. a proton and its adjacent electrons — has a linear dimension not 
larger than 

^Smax = C ■ Tact < 0.3 A 

(c: velocity of light), since the neutron-nucleus scattering dynamics during Tact 
cannot be causally influenced by other particles being more than Asmai apart 
from the nucleus. Consequently, the scattering system must necessarily be of "mi- 
croscopic" dimensions; it contains the scattering nucleus and a part of the adjacent 
electron density. 

However, this is not in conflict with the above estimate, for the following rea- 
son. As standard theory shows 1151 . S{q,uj) is peaked around the nuclear recoil 
energy Eq = t^q^ I2ra. The scattering time t^c is also given by the inverse of 
the width A_E of S{(i,uj), and S(iq,^) plays the role of the probability density 
distribution for transferring energy fioj from the neutron to the proton, when the 
momentum transfer is hq. That is, t^c ~ ti/ /S.E. (Interestingly, as Gidopou- 
los 1161 showed, t^c is also about the inverse of the energy spread of the pro- 
ton wavefunction after collision.) For a typical value AS ^ 10 eV, one gets 
Tsc ~ 10~^® — 10~^^ s. In other words, the scattering time Tsc gives a statistical 
measure of the length of the time interval during which an actual neutron-proton 
collision may occur — in the same way that the spatial extent of a particle wave- 
function gives a statistical measure of the extent of the region in which the particle 
may be found. 

To shed more light upon the issue of "relevant scattering time", one may also 
refer to the celebrated Margolus-Levitin theorem 1 17 1. Let us consider the neutron- 
proton system during the collisional process. Obviously, the initial and final states 
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of it are very different and so tliey can safely be assumed to be orthogonal. This 
theorem asserts that in takes at least a time Tx > (tt/i) / {2Es) for the system to 
evolve from its initial to any orthogonal final state. Es is the system's average en- 
ergy minus its ground state energy. T± provides a strict bound for the considered 
dynamical process 1 17,| . Note that in NCS one has Es ~ huj with huj taken at the 
peak center, and thus Es is of similar order as the aforementioned energy spread 
AE. Thus it is revealing that also this time T± is very similar to the aforemen- 
tioned scattering time, i.e. T± < Tsc- 



3 Theoretical model — NCS from open quan- 
tum systems 

In the following we present an outline of a recently proposed theoretical interpre- 
tation of the considered effect, which is based on the general theory of scattering 
from open quantum systems 1 18|. As explained above, the scattering system must 
necessarily be of "microscopic" dimensions (i.e. it is of the order of one A or less), 
and since it is embedded in condensed matter it represents an "open" quantum sys- 
tem 1211 . This point is crucial, since standard neutron scattering theory always 
assumes a condensed matter scattering system to be closed; see e.g. |T9ll20l . 

3.1 Scattering from closed systems 

First let us consider neutron scattering from a closed system consisting of N par- 
ticles with the same scattering length b, and the N-body Hamiltonian Htotai = 
Hq -\- V with the interaction 

V{v) = An(r), with A = h . (4) 

m 

m is the neutron mass, n(r) is the particle density operator 

1 

n{v)=yY.^{v^n,), (5) 

where V is the volume, and Rj is the position of the j-th particle; cf. the textbook 
(20). 

In the interaction picture, the Schrodinger equation is now (setting for simplic- 
ity h = 1) idt'^ — A n(r,t)4', with the perturbative solution 

= *(0) -iX [ n(r, t') dt'^{0). (6) 
Jo 

We write the transition probability W{t) between initial states tjji (with proba- 
bility Pi) and final states tpf of the scattering system to be given by 

Wit) = ^ I (^/ I A /' n(r, t') dt' I i,,) f p. (7) 

It should be noted that xpi and f are eigenstates of the A^-body Hamiltonian Hq 
omitting the probe system I19II20I . The transition probability is then given in the 



4 



form ^ ^ 

W{t) = [ dt' [ dt"Y,{^f I n{r,t') pn{v,t") \ ^f), (8) 
Jo Jo J 

with p = 5]]^ \ Tpi) Piii'i |, where we have noted that (r, t) = n(r, t). 

In an actual scattering experiment from condensed matter, we do not measure 
the cross-section for a process in which the scattering system goes from a specific 
initial state ipi to another state ipf, both being unobserved states of the many-body 
system. Therefore, one takes an appropriate average over all these states 11911201 . 
as done in Eq. l[7j. 

Given the initial (ko) and final (ki) momenta of an impinging neutron and 
introducing the momentum transfer q = ko — ki from the probe particle to the 
scattering system, the Fourier transform of the particle density reads 

n{r, t) = J rfqn(q, t) exp (i q • r) , (9) 

where, in the case of neutron scattering, cf. Eq. lO, 

n(q, i) = Yl '^^P [-'^ ■ • (10) 

Since n{r,t) is Hermitian, we have n^(q, t) — n(— q, t) and one obtains from 
Eq. 

W{t)^X^ f dt' f dt"Y,{^f\n{ci,t') pn{-^X)\^s)- (11) 
Jo Jo J 

Assuming as usual that '£f\tpf}{tpf \ = 1 [ !T9ll20l we get 

J2{i'fH(l,t')pn{-ci,t")\^f) =Tr[n{q,t')pn{~(i,t")] , (12) 
/ 

If the integration in Eq. i ll It is extended over all times (i.e., t oo), this 
ensues over-all energy conservation. This reproduces the well known result of 
standard neutron scattering theory, cf. 1191 1201 . Here, however, it is important 
to retain the finite duration of the scattering time, t < Tsc- This introduces an 
additional freedom into the theory, because we may be able to observe the influence 
of the decoherence on the scattering yield; see below. The result will be expressed 
in terms of the correlation function 

C(q,r) = rr[n(q,t)pn(-q,t + r)] = rr[n(q,0)pn(-q,r)] , (13) 

where we have utilized the fact that the scattering system is stationary. By intro- 
ducing the so-called scattering time Tsc, representing the time interval in which the 
scattering process may happen, we find 



W{rsc) = J dt' J dt"Ciq,t"-t') = J dt' J dri [C(g,r?) + C{q, -r])] . 



(14) 
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Here we use the stationarity of the correlation function 1201 . If we assume this 
function to be real, and that C{q, ??) ~ for > Tsc, we obtain the result 

W{tsc) ~2X\sc J dr,Ciq,Tj). (15) 



Now we can introduce the transition rate, W say, which is defined as 




Here the coiTelation function is analogous to the so-called intermediate function of 
neutron scattering theory |20|. This result for the scattering yield is analogous to 
that of standard theory. 



3.2 Dynamics of open systems and scattering 

We now consider the scattering system to be open and strongly interacting with 
its environment. We introduce a set of preferred coordinates { \^}}, cf. 1211 1221 
1231 , representing the relevant system's degrees of freedom coupled to the neutron 
probe. The density matrix p in ( |13t is then the reduced one in the space spanned 
by these states, and it is obtained by tracing out the degrees of freedom of the 
"environment" I21ll22l|23l . For simplify, throughout this section, let us denote 
this reduced density matrix by p too. 

In the subspace spanned by the preferred coordinates (also termed "pointer 
basis"), let us assimie a Lindblad-type equation | 24|(a) of the form 

dtp=-i[H,p] + Kp = £.p (17) 

which has the formal solution p(t) = e^*p(0). Let us look at a time-dependent 
expectation value 

{A{t)) = Tr [p{t)A) = Tr (e^V(O)A) = Tr (p(0)e^''A) , (18) 

where £t is defined by Tr {{CX) Y) = Tr {X {C^Y)). Thus it holds dtA{t) = 
A{t); cf. l24K a). Here is assumed that C is time independent. 
For the correlation functions like the one in Eq. l |13l ) it then holds 

{A{t)B) = Tr [p(0) (e^'*^) b] = Tr [Ae^' {Bp{0))] = Tr {ApB{t)) , 

(19) 

where psit), as defined in Eqs. ( 1191 ). obeys the equation dtpBit) = CpBit) with 
the initial condition ps(0) = Bp{0). Thus, except for the initial condition, we 
have to solve the same equation of motion as for the density matrix, Eq. dlVt . 

For simplicity of the further derivations, let us assume here a simple Lindblad- 
type ansatz for the master equation having only one Lindblad variable X. (In the 
real system we would have a multitude of such variables.) Thus 

dtp^^i[H,p]~K[X,[X,p]]^Cp, (20) 

where the constant K is real and K > 0, H is the reduced (or relevant) Hamilto- 
nian of a microscopic scattering system, and the double commutator term describes 
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decoherence (and/or dephasing). For simplicity of the further calculations, we fur- 
ther assume that we can take the preferred coordinates to commute with the total 
Hamiltonian 

H\o = £i\0, ^ U> = C U> • (21) 

This time evolution is now introduced into the correlation function C(q, r), 
Eq. il3l . A short straightforward calculation (see Ref. 1181 for full details) yields 
for the transition rate the result 

W = n J2 [-« (^«' - ^e) ^] exp [-K (C' - 0^ r] x 
Jo 

X (e ! n(-q, 0) I C'> I n(q, 0) p(0) U> dr . (22) 

The decoherence-free limit of this result (i.e. K = 0) corresponds to the con- 
ventional result of scattering theory. The oscillatory terms exp [—i {£^i — £()t] 
are due to the unitary dynamics caused by the commutator part —i [H, p] of the 
master equation l |20| l for p. These factors have the absolute value 1. If decoher- 
ence is present {K > 0), and especially if K^^ ~ Tsc, the additional contractive 
factors exp(-A: {£,' - (f r) < 1 can be seen to cause a decrease of the transition 
rate and thus of the associated cross-section. This can be illustrated as follows. 

Let us first assume that the reduced density operator p(0) can be chosen to 
be diagonal in the preferred ^—representation (which corresponds to the usual 
random phase approximation). Then each term of Eq. l l22t contains the factor 

(e|n(-q,0)|e')(eV(q,0)p(0)|O = KC I 'i(-q,0) I e'>|'(?|p(o)|0 > 0. 

(23) 

In the more general case with p{Q) being not diagonal in the ^—representation, one 
may proceed as follows. The decoherence factors exp(— if H' — t) imply that 
only terms with ^ ~ ^' contribute significantly to the transition rate. Thus we may 
conclude that, by continuity, all associated terms with ^ « ^' in Eq. l l22t should be 
positive, too. The further terms with ^ being much different from ^' can be positive 
or negative. But they may be approximately neglected, since they decay very fast 
and thus contribute less significantly to W\ cf. |22|. 

The main conclusion from the preceding considerations is that the time average 
over Tsc in Eq. l l22t always decreases the value ofW = W{tsc) /tsc, due to the 
presence of the contractive factors exp(— A' (5' — r) < 1. In other words, the 
effect of decoherence (and/or dephasing) during Tsc plays a crucial role and may 
lead to an "anomalous" decrease of the transition rate and the associated scatter- 
ing intensity. This result is in line with that of Ref. [25] , which investigated the 
standard expression of the double differential cross-section of neutron Compton 
scattering theory I14II15I by assuming decoherence of final and initial states of the 
scattering system. 

For further illustration let us consider the following two specific limiting cases: 
(A) For "vanishing" decoherence, if — > 0, the above contractive factors go to 1 
and thus the anomalous scattering effect disappears. I.e. the preceding result l l22t 
agrees with the conventional theoretical results |15||20| . (B) In the opposite case, 
K — > oo, only the "diagonal" terms with = ^' survive in Eq. l l22t and the related 
contractive factors go to 1; additionally, the oscillating factors become equal to 
1 . Consequently, the time-integration in Eq. M2\ has no effect and the rhs of this 
equation goes over to the standard expression Eq. l ll6b . Also this result is in line 
with conventional expectations. 
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4 Conclusions 



Theoretical considerations suggest the presence of attosecond quantum entangle- 
ment of the scattering protons and the surrounding electrons fT] [3] 1121 1131 1251 . 
Furthermore the usual Born-Oppenheimer approximation is not applicable 1131 
|16li25 26 1 . Moreover, recent NCS results from liquid HD and the equimolar H2- 
D2 mixture showing identical anomalous scattering, provided strong evidence that 
quantum exchange correlations between proton pairs cannot be the main physical 
reason for the effect {27]. For further theoretical discussions, see f25]. 

In contrast to standard scattering theory 1 19 20 1 of thermal and/or cold neu- 
trons, in which decoherence plays no role at all, our theoretical treatment of NCS 
is based on the physical fact that micro- and/or mesoscopic scattering systems are 
open quantum systems. This was shown to follow from the ultrashort scattering 
time of NCS and ECS. The revealed "anomalous" effect, which has no conven- 
tional interpretation, indicates that attosecond entanglement (and its decoherence) 
involving protons are quantum phenomena of broader significance and relevance 
than realized so far. 
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